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The origin of pseudospin symmetry and its broken in real 
nuclei are discussed in the relativistic mean field theory. In 
the exact pseudospin symmetry, even the usual intruder orbits 
have degenerate partners. In real nuclei, pseudospin symme- 
try is approximate, and the partners of the usual intruder 
orbits will disappear. The difference is mainly due to the 
pseudo spin-orbit potential and the transition between them 
is discussed in details. The contribution of pseudospin-orbit 
potential for intruder orbits is quite large, compared with 
that for pseudospin doublets. The disappearance of the pseu- 
dospin partner for the intruder orbit can be understood from 
the properties of its wave function. 
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Pseudospin symmetry was originally observed in spher- 
ical atomic nuclei more than 30 years ago , and later 
proved to be a good approximation in deformed nuclei 
including triaxial deformed nuclei This symme- 
try has been used to explain features of deformed nu- 
clei, including superdeformation and identical bands 
as well. The concept of pseudospin is based on ex- 
perimental observations that single particle states with 
non-relativistic quantum numbers (n,l,j = I + 1/2) 
and (n — 1,Z + 2,j = I + 3/2) lie very close in energy 
and can therefore be relabelled as pseudospin doublets: 
(n = n— 1, 1 = £ + 1, j = £±1/2). For example, the pair 
of orbits (3s 1 / 2 , 2d 3 / 2 ) can be viewed as a pseudospin 
doublet (2p!/ 2 , 2p 3 / 2 ). While, all pi/ 2 states are pseu- 
dospin singlets S1/2, and states with n = I, j = £ + 1/2 
are intruder orbital states without pseudospin partner. 

Since the suggestion of pseudospin symmetry in atomic 
nuclei, there have been comprehensive efforts to under- 
stand its origin. Apart from the rather formal relabeling 
of quantum numbers, various proposals for an explicit 
transformation from the normal scheme to the pseu- 
dospin scheme have been made in the last twenty years, 
and the pseudospin symmetry is proved to be connected 
with the special ratio in the strengths of the spin-orbit 
and orbit-orbit interactions J^-p^l- Despite the long his- 
tory of pseudospin symmetry, the origin of symmetry has 
eluded explanation. 

The relation between pseudospin symmetry and the 



relativistic mean field (RMF) theory Jl| was first noted 
in jl0|, in which Bahri et al. found that the RMF theory 
explains approximately the special ratio of the strengths 
of spin-orbit and orbit-orbit interactions in the non- 
relativistic calculations for exact pseudospin symmetry. 
No great progress was achieved until Ginocchio revealed 
that the pseudo-orbital angular momentum is nothing 
but the "orbital angular momentum" of the lower compo- 
nent of the Dirac spinor and built the connection between 
pseudospin symmetry and the equality in magnitude but 
difference in sign of the scalar potential V s (r) and vector 
potential V v (r) 14 Based on Relativistic Contin- 
uum Hartree-Bogoliubov (RCHB) [|l6],[l7) , it is shown 
that pseudospin symmetry is exact under a more general 
condition, d(V s + V v )/dr — 0, and the quality of the pseu- 
dospin approximation in real nuclei is connected with the 
competition between the pseudo-centrifugal barrier and 
the pseudospin-orbital potential |l^Jl^| . 

Here, in this paper, we will show why the pseudospin 
symmetry is approximate and broken in real nuclei, and 
explain the partner of intruder orbit is missing. We will 
first solve Dirac equation for two kinds of harmonic os- 
cillator potentials with exact spin symmetry and pseu- 
dospin symmetry, respectively. Taking 208 Pb as an ex- 
ample, we will investigate how the pseudospin symmetry 
is a good approximation for normal orbits and why in- 
truder orbit fails to have partner in real nuclei. A short 
summary is given at the end. 

The Dirac equation of a nucleon with mass M moving 
in an attractive scalar potential V s (r) and a repulsive 
vector potential V v (r) can be written as, 



■p + /3(M + V S ) + V V ] 



(1) 



For spherical nuclei, the nucleon angular momentum «7, 
and k = —fi(p ■ L + 1) commute with the Dirac Hamil- 
tonian, where (3, a, and L are respectively the Dirac 
matrix, Pauli matrix, and orbital angular momentum. 
The eigenvalues of k are n — ±(j + 1/2) with — for 
aligned spin (si/2, P3/2? etc.) and + for unaligned spin 
(Pi/2, d 3/2 , etc.). For a given n = ±1,±2, j = \k\ - ±, 

£ = \k+1/2\-~ ,1= |/t-l/2|-l/2. The wave functions 

can be classified according to their angular momentum 
j, k, and the radial quantum number n. In terms of the 
upper and lower radial functions F nK (r) and G nK (r) , Eq. 
(Ill), becomes M: 
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(l + ^)F nK (r) = (M- 



E nK -A) G nK (r), 



(2) 
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(— - -) G nK (r) = (M - E nK + V) F nK (r), (3) 
dr r 

with A = V v —V s and V = V s +V v . Eliminating G nK (r) or 
-PtikM one can get the following Schrodinger-like equa- 
tions, 



dr 2 



£{£ + !) 



{M + E nK -A)(M-E nK + V) 



+ 



dA d ^ 

dr dr r 



M + E n 



A 



F nK (r) = 0, (4) 



dr 2 



til + 1) 



-(M + E nK -A)(M-E nK +V) 



dV ^d k . 

dr dr r 
M — E nK + V 



G„«(r)=0, (5) 



and the same result can be obtained by solving either 
Eq.(§ or Eq.(§). 

For A = or — = 0, Eq.(|) is reduced to 



d 2 _ 1(1+1) 
dr 2 r 2 



F nK (r) 



= [(M + E nK - A)(M - E nK + V)] F nK {r), 



(6) 



The eigen energies, E nti , depend only on n and £, i.e., 
£; niK = E{n,£{£ + 1)) . For £ ^ 0, the states with j = 
£±1/2 are degenerate. This is a SU(2) spin symmetry. 

For the Dirac equation with the following vector and 
scalar potentials, 



V v (r) = V s (r) = X -Mu 2 r 2 , 



(7) 



A = V v (r) — V s (r) = 0, Eq.(||) can be solved analytically 
with the eigen energy, 



where Ai is a normalization constant and L^L 1 / 2 is the 
Laguerre polynomials. The lower radial wave function 
can be obtained from Eq.(||). For the spin partners, 
although the upper wave functions are the same, their 
lower wave functions are different because their k are 
different. 

In Fig. |l|, the energy spectrum for the Dirac equation 
with V v = V s = 2.5 r 2 is given by solving Eq.(0). Of 
dA 

course, A = or — — = means the missing of the spin- 
dr 

orbital interaction, which cannot be true for real nuclei. 

For V = or — — = 0, Eq.( p3) is reduced as 
dr 



dr 2 



£{£ + !) 



Gnuir) 



= [ (M + E nK - A)(M - E nK + V) } G nK (r). (10) 

The solution of the Dirac equation can be obtained from 
Eq.(|Io|) and its eigen energies, E nK , depend only on n and 
£, i.e., E nK = E(n, £(£ + 1) ). For I ^ 0, the states with 
j = £± 1/2 are degenerate. This is the pseudospin SU(2) 
symmetry. Except pseudospin singlet states with £ — 0, 
every aligned state ( j = £ + 1/2 = £ — 1/2) will have a 
degenerate unaligned partner ( j = £ — 1/2 = 1+1/2 ), 
and vice versa. Even the usual intruder orbital states 
may have partner now. These characters are very similar 
to the spin SU(2) symmetry. 

For the Dirac equation with the following vector and 
scalar potentials, 



V v {r) = -V s (r) = X -Muj 2 r 2 , 



(11) 



V = V v (r) + V s (r) = 0, Eq.(|g) can be solved analytically 
with the eigen energy, 



u;(2n + £~-) (n=l,2,...). 



(12) 



The states with the same n and £ will be degenerate. The 
corresponding lower radial wavefunction is, 



cu{2n + £--) (n=l,2,...). 



(8) 



Now states with the same n and £ will be degenerate. Of 
course as the potential is a special harmonic oscillator, 
the spectrum is more degenerated. The corresponding 
upper radial wavefunction is, 



F nK (r) = Ai (ar) 



1 2 2 

!+i e ~2 a T L l+1 / 2 t~. 2 ~ 2 \ 



n-l 



>V), 



( a 2 = ui\J M(E nK - M) ), 



(9) 
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G nK (r) = A 2 (ar) e+1 e~ 2~ ' L n J(a 2 r 2 ) 

( a 2 = uy/M(E nK + M) ), (13) 

where A2 is the normalization constant and L^ 1 / 2 is the 
Laguerre polynomials. The upper radial wave function 
can be obtained from Eq.(||). For the pseudospin part- 
ners, although the lower wave functions are the same, 
their upper wave functions are different due to k. 
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In Fig.0, the energy spectrum for the Dirac equation 



with V v = —V s = 2.5 r is given by solving Eq.(|10[). 
Similar to Fig.jl], the spectrum is more degenerate due 
to the special harmonic oscillator used here. To make 
sure that the radial quantum number in Figj^ is correct 
and even the intruder orbits have pseudospin partners, 
the spectrum in Figj2| is confirmed by solving the Dirac 
equation ([[]) numerically as in Refs. |I6| , |l7| . 

Now the questions arise: in real nuclei why the pseu- 
dospin symmetry is broken and there is no partner for 
intruder orbit. 

In real nuclei, — — 7^ 0, and pseudospin symmetry 
dr 

is only an approximation. The quality of the pseu- 
dospin symmetry approximation depends on the compe- 
tition between the contributions of the pseudo-centrifugal 
1(1 + 1) 

potential, ^ — , and the pseudospin-orbit potential, 

dV re „ i-k 

in Eq. (|). 



dr r(M — E + V)' 

In [ ^8|Jl9| , Meng et al. analyzed the competition of 
the two potentials for pseudospin partners, and found 
that the contribution of the pseudospin-orbit potential is 
small, compared with that of the pseudo-centrifugal po- 
tential. Following the arguments and conjectures in Ref. 
IP^P^I , pseudospin symmetry is discussed for the nuclei 
with deformed shape or isospin asymmetry 21 23). 

To understand why the intruder orbit has no partner 
is a challenging problem. Lots of works have been done 
along this line. In Ref. |^J,^5|, the wave functions of 
the pseudospin partners and the structure of the radial 
nodes have been investigated. Here combining the energy 
spectrum in Fig]| with that obtained from RCHB [ |l6| , |l7| , 
we will try to understand why the intruder orbit has no 
partner. 

In Fig|s, multiplied with the wave function G 2 , the ef- 

re dV 

fective pseudospin-orbit potential — — and the effective 



dr 



1) 



in 208 Pb have 



centrifugal barrier (M — E + V)- 

7— 

been compared in detail for the pseudospin singlets, the 
pseudospin doublets, and the intruder orbits. 

For pseudospin singlets, we choose lpi/2 and 2p 1 / 2 as 
examples, and find that the pseudo-centrifugal potential 
is always zero. For the usual pseudospin doublets, the 
contribution of the effective pseudospin-orbit potential is 
much smaller than that of the effective pseudo-centrifugal 
potential, and the pseudospin symmetry is a good ap- 
proximation. For the intruder orbits, the contribution 
of the effective pseudospin-orbit potential is comparable 
with and sometimes even larger than that of the effec- 
tive pseudo-centrifugal potential. It can be seen that the 
pseudospin-orbit potential is responsible for the broken 
of pseudospin symmetry. So far it is not explained where 
the pseudospin partner of the intruder orbit in real nuclei 
has gone. In the following we will try to get the answer 
from the properties of the wave function. 

For bound states, their solution can be obtained from 



either one of the Schrodinger-like equations, Eqs.(Q) and 
(|). Let F(r) = r~ K f(r), from Eq.(|) one get: 



^- = r*(M + E nK -A) G. 



(14) 



For r ~ 0, the solution of bound states for Eqs.(|2]) and 
(0) can be obtained analytically as: 



F(r) = A Sign(re) X Hz + M ^ J\ K+ i/ 2 \(kr) (15) 



2M + V - A 



G(r) = A 



E nK -M-V 
2M + V-A 



VrJ\ K ~i/2\(kr), 



(16) 



where k = y/-(M + E- A) (A/ -E + V) and A is a nor- 
malization factor ( here assumed to be positive ). For 
re > 0, if the node number of lower wave function G is 
zero: tiq — 0, then G will be always positive. As r ~ 0, 
one has: 



Mr) =r n (- 



+ OH 2 ), 



T(l+n) 

and depending on re, f(r) can be written as: 



(17) 



f(r) ~ A 



E„ 



M 



2 -|k+1/2| 



2M + V-A T(l 



1/21) 



r 2K+1 re > 
-1 re < 



(18) 



As the factor M + E nK — A in Eq.(n4) is always posi- 

df 

tive, for no = and re > 0, — will be always positive, 

dr 

i.e., f(r) increases with r. While for bound states, the 
wave function at large r should decrease as an exponen- 
tial function of r, i.e., /(oo) = 0. Therefore for rig = 
and re > 0, the boundary condition of / for bound state 
can not be satisfied and there is no corresponding bound 
states. This means that the pseudospin partner of the in- 
truder orbit will disappear. Similar conclusion has been 
drawn in Ref. [ p5|j24| . To restore the pseudospin symme- 
try for intruder orbit, one possibility is that the factor 
(Af + E nK — A) can change sign with r, as the case in 
Fig-1 

In summary, the origin of pseudospin symmetry and its 
broken in real nuclei are discussed in the RMF theory. In 
the exact pseudospin symmetry, all the orbits except the 
singlets P1/2 have partners. In real nuclei, pseudospin 
symmetry is approximate, and the partners of the usual 
intruder orbits will disappear. The competition between 
the pseudo-centrifugal and pseudospin-orbit potentials 
decides the quality of pseudospin symmetry. The contri- 
bution of pseudospin-orbit potential for intruder orbits is 
quite large, compared with that for pseudospin doublets. 
The disappearance of the pseudospin partner for the in- 
truder orbit can be understood from the properties of its 
wave function. 

Valuable discussions with Prof. J.Y. Zeng are grate- 
fully acknowledged. This work was partly supported by 



3 



the Major State Basic Research Development Program 
Under Contract Number G2000077407 and the National 
Natural Science Foundation of China under Grant No. 
10025522, 10047001, and 19935030 . 



FIG. 1. The eigen energy for the Dirac equation with 
V v (r) = V 3 {r) = ^MwV, M = 10.0, and u = 1.0. 
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FIG. 2. Same as Fig.§ but for V v (r) = -V s (r) = ^Muj 2 r 2 . 

(a) Pseudospin singlet states 

(b) Pseudospin doublet states 

(c) Intruder orbits 

FIG. 3. The competition between the effective 
pseudocentrifual G 2 (r)(E - M - V)£(£ + l)/r 2 ( solid line 

) and pseudospin-orbital potentials —G 2 (r)^- — ( dashed 

dr r 

line ). 
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